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1.  INTRODUCTION 

Until  recently  the  term  Interior 
Ballistics  referred  to  the  computation 
of  head  end  pressure  and  thrust  in 
solid  propellant  rocket  motors  by 
quasi-steady  analysis.  However,  with 
recently  developed  technology1*2,  it  is 
now  possible  to  predict,  given  the 
necessary  propellant  and  nozzle 
admittance  and  response  parameters, 
the  linear  stability  of  a rocket  motor 
subjected  to  small,  combustion 
amplified,  disturbances,  i.e., 
combustion  instability.  With  the 
advent  of  low  signature  constraints, 
this  capability  has  assumed 
considerable  importance  because  of 
the  high  probability  of  combustion 
instability  in  low  signature  motors 
and  the  attendant  jeopardy  to 
successful  motor  (and  system) 
operation  that  instability  creates. 

A major  shortcoming  of  current 
combustion  instability  prediction 
technology  is  that  it  is  computationally 
divorced  from  conventional  interior 
ballistic  analysis.  This  makes 
instability  analysis  difficult  because 
interior  ballistic  data  must  be  hand 
loaded  into  the  stability  codes. 
Moreover,  it  has  fostered  a functional 
split,  with  interior  ballistic  analysis 
the  domain  of  the  grain  designer  while 
combustion  instability  analysis  is  that 
of  the  "instability  specialist”.  Due  to 
the  difficulties  associated  with 


running  combustion  instability 
analyses  and  the  aforementioned 
functional  split,  combustion 
instability  has  not  been  an  integral 
part  of  the  motor  design  process. 
Clearly,  a successful  motor  design  for 
any  low  signature  application  should 
include  a combustion  instability 
prediction  aspect. 

The  remedy  for  this  shortcoming 
would  seem  to  be  obvious:  combine 
conventional  interior  ballistics 
analysis  and  stability  codes  into  a 
single  complete  interior  ballistics  code 
that  the  motor  designer  can  employ, 
i.e.,  a code  that  performs  pressure, 
thrust,  and  linear  stability  margin 
versus  time  analyses  virtually 
simultaneously.  In  this  way,  it  is  easy 
for  the  motor  designer  to  consider 
combustion  instability.  It  is  with  this 
concept  in  mind  that  the  development 
of  the  computer  program  described 
herein  was  undertaken.  Its 
development  is  based  upon  the 
utilization  of  the  previously  developed 
three-dimensional  finite  element 
combustion  instability  prediction 
code,  FLAP32.  The  utilization  of 
FLAP3  requires  automatic  finite  ele- 
ment mesh  and  boundary  condition 
generation,  and  the  coupling  of  FLAP3 
with  an  interior  ballistics  analysis 
code,  therefore,  requires,  in  effect,  a 
coupling  of  the  ballistics  code  output 
with  the  mesh  generator  input.  This 
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process  is  the  essence  of  a development 
of  a combined  code. 

2.  THREE-DIMENSIONAL 
FINITE  ELEMENT  MESH 
GENERATION 

FLAP3  (Fluid  Analysis  Program,  3 
Dimensions)  performs  a linear 
acousto-modal  analysis  of  the 
irrotational.  motions  of  an  inviscid, 
compressible  fluid  coupled  to  the 
irrotational  motions  of  a nearly 
incompressible,  linearly  viscoelastic 
solid,  and  a linear  potential  flow 
anaylsis  of  the  irrotational  motions  of 
an  inviscid,  imcompressible  fluid,  and 
then  determines  the  effect  of  the  flow 
field  and  of  combustion  on  the 
calculated  acoustic  oscillations2--'.  This 
combustion  instability  analysis  is 
performed  at  different  points  in  time, 
beginning  at  initial  combustion,  or 
time  zero.  The  output  of  FLAP3  is 
modal  frequency  and  an  evaluation  of 
the  pressure  growth/decay  coefficient 
a for  each  mode  of  vibration  analyzed. 
A positive  net  value  of  a indicates  a 
growth  of  pressure  oscillations  and, 
therefore,  instability  while  a negative 
value  of  a is  an  indication  of  decaying 
oscillations,  or  stability. 

The  FLAP3  analysis  utilizes  the 
finite  element  method  and  models  the 
acoustic  cavity  and  propellant  grain 
as  an  assemblage  of  three-dimensional 
quasi-hexahedral  elements  connected 


at  the  corner  nodes.  The  mesh 
generation  code  used  in  conjunction 
with  FLAP3  is  FLESH3  (Fluid  Mesh,  3 
Dimensions)  which  was  developed  for 
the  purpose  of  generating  input  data 
for  FLAP3.  A detailed  description  of 
the  use  of  FLESH3  is  found  in 
Reference  3.  FLESH3  generates  the 
numbered  nodal  points  and  nodal 
coordinates  and  identifies  each  node 
as  to  whether  it  lies  in  the  acoustic 
cavity,  on  the  cavity-grain  interface,  or 
in  the  grain;  generates  the  quasi- 
hexahedral  elements,  designated  by 
the  eight  numbered  nodal  points 
defining  each  element,  and  identifies 
each  element  as  to  whether  it  is  a 
cavity  element,  a cavity  element 
adjacent  to  the  cavity-grain  interface, 
or  a solid  propellant  element;  and 
generates  the  cavity-grain  interface 
element  surfaces  (burning  surfaces) 
and  identifies  each  as  to  the  direction 
of  the  surface  normal.  This 
information,  along  with  a few 
additional  input  data  relative  to  gas 
and  propellant  properties  and  type  of 
acoustic  mode(s)  to  be  analyzed,  is 
necessary  and  sufficient  to  activate  a 
combustion  instability  analysis  by 
FLAP3,  given  adequate  computa- 
tional facilities. 

Input  to  FLESH3  is  in  the  form  of  a 
definition  of  global  curves,  defining 
cavity  and  grain  boundaries;  a 
designation  of  blocks  of  points, as  to 
whether  they  are  cavity,  interface  or 
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grain  points;  a designation  of  blocks  to 
be  generated  ,by  their  nodal  indices; 
and  a designation  of  the  number  and 
location  of  longitudonal  cross- 
sections,  defining  the  number  of  layers 
of  elements.  The  objective  is  that  of 
minimizing  the  amount  of  input  to 
FLESH3.  This  can  be  done  by  the  use 
of  an  executive  routine  which  accepts 
as  input  a small  number  of  geometrical 
parameters  from  the  interior  ballistics 
code  and,  in  turn,  activates  FLESH  3. 
The  development  of  this  routine  is  the 
major  thrust  of  this  report. 

3.  INTERIOR  BALLISTICS 
CODE 

The  interior  ballistics  analysis  code 
used  in  the  formulation  was  developed 
by  Aerojet  Solid  Propulsion 
Company4.  The  FLESH3  geometrical 
input  parameters  are,  therefore,  those 
presented  in  that  code,  but  similar 
parameters  would  be  obtained  from  a 
consideration  of  a comparable  interior 
ballistics  code.  The  description  of  the 
developed  executive  routine  will  be 
related  to  the  Aerojet  code,  but  the 
procedure  followed  in  the  development 
is  general. 

4.  DEVELOPMENT  OF  THE 
EXECUTIVE  ROUTINE 

A listing  of  the  formulated  computer 
code  GRNMSH  (Grain  Mesh)  is  found 
in  Appendix  A.  The  code  is  based  upon 


the  star  design  found  in  Reference  4 
and  reproduced  in  Figure  1.  One 
computer  card  of  input  parameters 
(seven  values)  automatically 
generates  the  three-dimensional  finite 
element  mesh  and  boundary 
conditions  needed  for  a FLAP3 
analysis.  Due  to  the  dihedral 
symmetry  provision  in  FLAP3,  only 
the  smallest  repeating  segment  need 
be  analyzed.  This  segment,  for  the 
general  star  geometry,  is  shown  in 
Figures  2 and  4;  with  the  necessary 
geometry  parameters  identified.  The 
only  difference  between  Figure  2 and 
Figure  4 is  that  the  angle  $ is  zero  in 
Figure  4.  This  condition  necessitated  a 
slightly  different  internal  provision  in 
GRNMSH.  Figures  3 and  5 indicate:  (a) 
the  key  global  curves  (others  were 
generated,  as  can  be  seen  from  the 
included  example),  and  (b)  the  key  part 
indices.  A working  knowledge  of 
FLESH33  is  necessary  in  order  to 
understand  the  intricacies  of 
GRNMSH,  and  to  modify  and  add  to  it, 
but  not  in  order  to  use  it. 

The  use  of  GRNMSH  can  best  be 
understood  by  following  an  example. 
Table  1,  from  Reference  4,  lists  the 
twelve  cases  of  the  star  design  upon 
which  the  formulation  of  the  executive 
routine  is  based.  The  generated 
meshes  for  the  twelve  cases  are  shown 
in  Figures  6 through  9.  The  special 
case  of  a shell  design  (circular 
cylinder),  for  which  the  parameter  <t>  is 
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Figure  1.  Star  design  with  outside  round  on  propellant  tip;  N=4 
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Figure  3.  Cross-section  of  repeating  segment  of  star  design  showing 
important  input  to  mesh  generator  (FLESH  3):  (a)  Global  curves 
and  (b)  Part  indices  (matches  Figure  2). 
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Figure  4.  Cross-section  of  repeating  segment  of  star  design;  independent 
parameters:  N,  Rc,  w,  Ya,  a , L;  4>  = 0 


Figure  5.  Cross-section  of  repeating  segment  of  star  design  showing 
important  input  to  mesh  generator  (FLESH3):  (a)  Global  curves 
and  (b)  Part  indices  (matches  Figure  4). 
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TABLE  1.  STAR  DESIGN  PARAMETERS-SYMMETRY  7 

(N  = 7) 


STAR 

(Design  No.) 

w 

In. 

a 

deg 

deg 

Y. 

In. 

Yc 

In. 

L 

in. 

D-5a 

10 

3.5 

80 

0.0 

El 

0.075 

D-5b 

10 

3.5 

85 

0.0 

mm 

0.068 

D-5c 

10 

3.5 

75 

0.0 

0.2 

0.136 

D-6a 

10 

3.5 

85 

0.0 

0.5 

0092 

4 1 

D-6b 

10 

3.5 

80 

0.0 

0.5 

0.118 

D-6c 

10 

3.5 

75 

0.0 

0 5 

0095 

D-7a 

10 

3.5 

85 

2.5 

02 

0 043 

4 8 

D-7b 

10 

3.5 

80 

2.5 

0.2 

0 084 

4 1 

D-7c 

10 

3.5 

75 

2.5 

0 2 

0047 

D-8a 

10 

3.5 

85 

2.5 

0.5 

0.092 

D-8b 

10 

3.5 

80 

2 5 

0.5 

0.035 

D-8c 

10 

3.5 

75 

2.5 

0.5 

0.069 

L3jl 

Flgur«6.  Cross-section  showing  generated  finite  element  mesh  for  star 
design  D-5  (see  Table  1). 


Figure  7.  Cross-section  showing  generated  finite  element  mesh  for  star 
design  D-6  (see  Table  1). 


Figures.  Cross-section  showing  generated  finite  element  mesh  for  star 
design  D-7  (see  Table  1). 
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Figure  10.  Cross-section  showing  generated  finite  element  mesh  for  shell 
design. 
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input  equal  to  ^VN,  as  generated  by 
GRNMSH  is  shown  in  Figure  10. 

Sample  computer  input  and  output 
for  the  star  design  geometry  D-5a 
( Table  1)  is  found  in  Appendix  B.  In 
sequential  order  is:  the  input  to 
GRNMSH,  the  output  from  GRNMSH 
(or  the  input  to  FLESH  3),  and  the 
output  from  FLESH3.  Again, 
familiarity  with  FLESH3,  by  way  of 
Reference  3,  will  enhance  one’s 
understanding  of  the  output  from 
FLESH3. 

It  can  be  noted  that  the  length  of  the 
grain,  and  provision  for  five 
equidistant  sections  were  internally 
set  in  the  program,  but  could, 
alternatively,  be  input  as  an 
additional  set  of  parameters.  It  can 
also  be  noted  that  the  value  of  e , the 
necessary  offset  radius  of  the  nodes 
closest  to  the  axis  origin,  is  internally 
set  at  0.1  inches. 


t 

5.  CONCLUSIONS 

A 

An  executive  routine  has  been 
developed  which  couples  an  interior 
ballistics  analysis  code  with  a 
combustion  instability  analysis  code; 
thus,  in  concept,  providing  the  rocket 
motor  designer  with  the  potential  for  a 
more  inclusive  and,  therefore,  a 
superior  design  approach. 


Basically,  the  proof  of  concept  has 
been  established  and  a fundamental 
executive  routine  developed  — one 
which  can  be  built  upon.  The  format  of 
the  program  is  simple  and  additional 
curves  defining  the  progressive  burn 
geometry  of  the  grain  can  be  easily 
added  to  the  routine.  The  development 
is  conceptually  complete;  the 
additional  effort  will  be  that  of 
enlarging  the  code  through  repetitive 
extension. 
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APPENDIX  A 

LISTING  OF  THE  COMPUTER  PROGRAM  GRNMSH 
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APPENDIX  B 

EXAMPLE  THREE-DIMENSIONAL  FINITE  ELEMENT  MESH 
GENERATION  FOR  STAR  DESIGN  (D-SA) 
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